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Eclipse Mapping: 

Astrotomography of Accretion Discs 
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Abstract The Eclipse Mapping Method is an indirect imaging technique that trans¬ 
forms the shape of the eclipse light curve into a map of the surface brightness dis¬ 
tribution of the occulted regions. Three decades of application of this technique to 
the investigation of the structure, the spectrum and the time evolution of accretion 
discs around white dwarfs in cataclysmic variables have enriched our understand¬ 
ing of these accretion devices with a wealth of details such as (but not limited to) 
moving heating/cooling waves during outbursts in dwarf novae, tidally-induced spi¬ 
ral shocks of emitting gas with sub-Keplerian velocities, elliptical precessing discs 
associated to superhumps, and measurements of the radial run of the disc viscosity 
through the mapping of the disc flickering sources. This chapter reviews the princi¬ 
ples of the method, discusses its performance, limitations, useful error propagation 
procedures, and highlights a selection of applications aimed at showing the possible 
scientihc problems that have been and may be addresses with it. 


1 Context and Motivations 

In weakly-magnetic cataclysmic variables (CVs) a late-type star overhlls its Roche 
lobe and feeds a companion white dwarf (WD) via an accretion disc, which usually 
dominates the ultraviolet and optical light of the system iol. The temperatures in 
CV discs may vary from 5000 K in the outer regions to over 50000 K close to the 
disc centre, and the surface density may vary by equally signihcant amounts over 
the disc surface. Therefore, the spectrum emitted by different regions of the accre¬ 
tion disc may be very distinct. Furthermore, aside of the component stars and the 
accretion disc, additional sources contribute to the integrated light from the binary: 
a bright spot (BS) or stream forms where the gas from the mass-donor star hits the 
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outer edge of the disc, tidally-induced spiral structures might develop in the outer 
regions of extended discs ESllll, and a fraction of the transferred mass may also be 
ejected from the binary in a wind from the disc surface ||9] . Because what one 

directly observes is the combination of the spectra emitted from these diverse re¬ 
gions and sources, the interpretation of CVs observations is usually plagued by the 
ambiguity associated with composite spectra. The most effective way to overcome 
these difficulties is with spatially resolved data, in which the light from the different 
sources might be disentangled. However, at typical sizes of less than a solar radius 
and distances of hundreds of parsecs, accretion discs in CVs are seen at angular di¬ 
ameters of tens of micro-arcseconds, and spatially resolving them is well beyond the 
reach of current direct imaging interferometric techniques, both at optical and radio 
wavelength ranges. Thus, resolving an accretion disc in a CV is presently possible 
only via indirect imaging. 

Developed in the 1980’s by Horne the Eclipse Mapping Method is an indi¬ 
rect imaging technique which provides spatially resolved observational constraints 
of accretion discs in CVs on angular scales of micro-arcseconds. It assembles the 
information contained in the shape of the eclipse into a map of the accretion disc 
surface brightness distribution. The following sections describe the technique and 
provide a set of examples aiming at illustrating the wealth of possible applications. 


2 Principles & Inner Workings 


The three basic assumptions of the standard eclipse mapping method are; (i) the 
surface of the secondary star is given by its Roche equipotential, (ii) the bright¬ 
ness distribution is constrained to the orbital plane, and (iii) the emitted radiation 
is independent of the orbital phase. While assumption (i) is reasonably robust and 
always employed, the others are simplifications that do not hold in all situations and 
may be relaxed with three-dimensional (3D) eclipse mapping implementations (see 
SectEl. 

A grid of intensities centred on the WD, the eclipse map, is defined in the or¬ 
bital plane. One usually adopts the distance from the disc centre to the internal 
Lagrangian point, Rlb as the length scale. With this definition the primary lobe has 
about the same size and form for any reasonable value of the binary mass ratio q 
{=M 2 /M\ , where M 2 and M\ are the masses of the mass-donor star and the WD, re¬ 
spectively) ||33]. If the eclipse map is an N points flat, square grid of side ARli, each 
of its surface element (pixel) has an area /N and an associated intensity /,. 

The solid angle comprised by each pixel as seen from the earth is then 


02 


(ARli)2 1 ' 
N if 


COS i , 


( 1 ) 


where d is the distance to the system. The value of X defines the area of the eclipse 
map while the choice of N sets its spatial resolution. 
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The model eclipse light curve m(0) is derived from the intensities in the eclipse 
map, 


m 


N 


( 0 ) = 02 ^ . 

j=l 


( 2 ) 


The eclipse geometry specifies the fractional visibility of each pixel as 

a function of orbital phase and may include fore-shortening and limb darkening 
factors llT5ll50ll65]l . The fractional visibility of a given pixel is obtained by dividing 
the pixel into smaller tiles and evaluating the Roche potential along the line of sight 
for each tile to see if the potential falls below the value of the equipotential that 
defines the Roche surface of the mass-donor star. If so, the tile is occulted. The 
fractional visibility of the pixel is then the sum of the visible tiles divided by the 
number of tiles. 

The eclipse geometry is determined by the inclination i, the binary mass ratio q 
and the phase of inferior conjunction 0o Ha ESI- These parameters set the shape 
and extension of the projected shadow of the mass-donor star in the orbital plane. 
As the binary phase (j) changes during eclipse, the resulting shadow rotates around 
the LI point, progressively covering/uncovering different parts of the accretion disc. 
This creates a grid of criss-crossed ingress/egress arcs in the orbital plane. A pixel 
with coordinates {x,y) within the region covered by this grid disappears and reap¬ 
pears from eclipse at a particular pair of binary phases , 0 Thus, the eclipse 

geometry sets the connection between the image space (the eclipse map) and the 
data space (the light curve). 

Figure[T] shows an example of eclipse geometry and depicts the connection be¬ 
tween the image and data spaces. Information about the location of the occulted 
brightness sources is embedded in the shape of the eclipse light curve. The width 
of the eclipse increases as the light source moves closer to the LI point, whereas 
it is displaced towards negative (positive) phases if the light source is in the lead¬ 
ing (trailing) disc side. For example, the phase width and range of the eclipse in 
the upper right panel of Fig.[T]tells us that the corresponding light source (the blue 
pixel in the left-hand panel of Fig.[TJ is located in the leading, near side of the disc. 
Figured illustrates the simulation of the eclipse of a fitted brightness distribution 
while showing the comparison between the resulting model light curve and the data 
light curve. The geometry in this case is i = 81° and q = 0.5. The left-hand pan¬ 
els show the data light curve (small dots) and the model light curve (solid line) as 
it is being drawn at five different orbital phases along the eclipse (indicated in the 
upper right corner). The right-hand panels show the best-fit disc brightness distribu¬ 
tion and how it is progressively covered by the dark shadow of the mass-donor star 
during the eclipse. The slope of the light curve steepens whenever a bright source 
is occulted or reappears from eclipse. The two asymmetric arcs of the brightness 
distribution of this example are occulted at distinct phases, leading to a V-shaped 


* Pixels outside the region covered by the grid of criss-crossed arcs are never eclipsed, and have 
no corresponding ingress/egress phases (p,',pj). Accordingly, there is no information about the 
surface brightness distribution of these unocculted regions in the shape of the eclipse. 
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Fig. 1 Left-hand panel: A grid of criss-crossed ingress/egress arcs for the eclipse geometry / = 81°, 
q = 0.15, and an eclipse map with X = 2. The ingress/egress arcs for the blue and red pixels are 
shown in the respective colour. A dotted line depicts the primary Roche lobe; a cross marks the disc 
centre. The mass-donor star is to the right of the panel and the stars rotate counter-clockwise (equiv¬ 
alently, the observer and the mass-donor star shadow rotate clockwise). Middle-panels: grayscale 
brightness distributions with a faint, extended disc plus bright, narrow Gaussian spots at the posi¬ 
tions of the blue and red pixels in the left panel. Right-hand panels: the model eclipse light curves 
obtained by convolving the eclipse geometry in the left panel with the brightness distributions in 
the middle panels. Vertical dotted lines mark the ingress/egress of the disc centre; blue/red vertical 
dashed lines depict the ingress/egress phases of the blue/red pixels. 


eclipse with two bulges at ingress and at egress. The flux at phase 0=0 does not go 
to zero because a sizeable fraction of the disc remains visible at mid-eclipse. 

Given the eclipse geometry, a model light curve can be calculated for any as¬ 
sumed brightness distribution in the eclipse map. A computer code then iteratively 
adjusts the intensities in the map (treated as independent parameters) to And the 
brightness distribution the model light curve of which fits the data eclipse light curve 
within the uncertainties (Fig.O. Because the one-dimensional data light curve can¬ 
not fully constrain a two-dimensional map, additional freedom remains to optimise 
some map property. A maximum entropy (MEM) procedure ll56l l57l is used to se¬ 
lect, among all possible solutions, the one that maximises the entropy of the eclipse 
map with respect to a smooth default map. 

The entropy of the eclipse map p with respect to the default map q is defined as 
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Fig. 2 Simulation of a disc eclipse (i = 81°, ^ = 0.5). Left-hand panels: data light curve (dots) and 
model light curve (solid line) at five different binary phases (indicated in the upper right comer). 
Vertical dotted lines mark mid-eclipse and the ingress/egress phases of the disc centre. Right-hand 
panels: visible parts of the best-fit eclipse map (in a false colour blackbody logarithmic scale) at the 
corresponding binary phases. Dashed lines show the primary Roche lobe and the ballistic trajectory 
of the gas from the mass-donor star; crosses mark the centre of the disc. The mass-donor star is to 
the right of each panel; the stars and the accretion disc gas rotate counter-clockwise. 
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(3) 


where p and q are written as 


Pi = 


Ij 

Li: 4 


^j = 


Dj 

LkDk' 


(4) 


Other functional forms for the entropy appear in the literature lll5l l60ll . 
equivalent to (12l when p and q are written in terms of proportions. 

The default map Dj is generally defined as a weighted average of the 
in the eclipse map, 

_ Lk Oijkh 
■’ Li (Ojk ’ 


These are 
intensities 

(5) 


where the weight function cOjk is specified by the user. A priori information about 
the disc (e.g., axi-symmetry) is included in the default map via (Ojic- Prescriptions 
for the weight function (Ojk and their effect in the reconstructions are discussed in 
Sect. [3 In the absence of any constraints on Ij, the entropy has a maximum S„,ax = 0 
when pj = qj, or when the eclipse map and the default map are identical. 

The consistency of an eclipse map may be checked using the a constraint 
function, 


X 


2 


]_ y [ 

M V / 


L K0) 
0= 


( 6 ) 


where c/(0) is the data light curve, (7(0) are the corresponding uncertainties, r(0) 
is the residual at the orbital phase 0, and M is the number of data points in the light 
curve. Alternatively, the constraint function may be a combination of the x^ and the 
R-statistics ca, 

1 M-l 

R=^== ^ r(0)r(0 + l), (7) 

VM-1 

to minimise the presence of correlated residuals in the model light curve m. For 
the case of uncorrelated normally distributed residuals, the R-statistics has a Gaus¬ 
sian probability distribution function with average zero and unity standard deviation. 
Requiring the code to achieve an R = 0, is equivalent to asking for a solution with 
uncorrelated residuals in the model light curve. 

The final MEM solution is the eclipse map that is as close as possible to its de¬ 
fault map as allowed by the constraint imposed by the light curve and its associated 
uncertainties llTSl l50l . In mathematical terms, the problem is one of constrained 
maximisation, where the function to maximise is the entropy and the constraint is a 
consistency statistics that measures the quality of the fitted model to the data light 
curve. Different codes exist to solve this problem. Many of the eclipse mapping 
codes are based on the commercial optimisation package MEMSYS Alterna¬ 
tive implementations using conjugate-gradients algorithms lfT4l [la, CLEAN-like 
algorithms ll58l and genetic algorithms are also used. 
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Rutten et al. El found that the entropy function can be a useful tool to signal 
and to isolate the fraction of the total light which is not coming from the accretion 
disc plane. They noted that when the light curve is contaminated by the presence 
of additional light (e.g., from the mass-donor star) the reconstructed map shows 
a spurious structure in its outermost region. This is because the eclipse mapping 
method assumes that all the light is coming from the accretion disc, in which case 
the eclipse depth and width are correlated in the sense that a steeper shape corre¬ 
sponds to a deeper eclipse. The addition of an uneclipsed component in the light 
curve (i.e., light from a source other than the accretion disc) ruins this correlation. 
To account for the extra amount of light at mid-eclipse and to preserve the bright¬ 
ness distribution derived from the eclipse shape the algorithm inserts the additional 
light in the region of the map which is least affected by the eclipse. Since the entropy 
measures the amount of structure in the map, the presence of these spurious struc¬ 
tures is flagged with lower entropy values. The correct offset level may be found by 
comparing a set of maps obtained with different offsets and selecting the one with 
highest entropy. Alternatively, the value of the zero-intensity level can be included 
in the mapping algorithm as an additional free parameter to be fitted along with the 
intensity map in the search for the MEM solution |l9]|49l. A detailed discussion on 
the reliability and consistency of the estimation of the uneclipsed component can be 
found in iflhl . 

Figure[3]gives an example of the convergence process of a MEM reconstruction. 
The code starts from a flat map and quickly evolves towards an axi-symmetric Gaus¬ 
sian map which reproduces the gross features of the data light curve. However, at 
this point the model light curve is not a good match to the data = 16.7), failing 
to reproduce the double stepped ingress/egress bulges of the eclipse shape. Several 
more iterations are required in order to match these features in the data light curve, 
which demands building two asymmetric bright arcs in the leading and trailing sides 
of the disc. The reconstruction converges for a final unity value after 233 itera¬ 
tions. 


3 Performance & Limitations 

A crucial aspects of eclipse mapping is the selection of the weight function for the 
default map, cojk, which allows the investigator to steer the MEM solution towards 
a determined type of disc map. The choice (Ojk = 1 results in a uniform default 
map which leads to the most uniform eclipse map consistent with the data. This 
happens not to be a good choice for eclipse mapping, because it results in a map 
severely distorted by criss-crossed artefacts ifT^I^ISSll as the flux of point sources 
is spread along their ingress and egress arcs (Eig.|4ja)). This led to the adoption of 
a weight function which sets the default map as an axi-symmetric average of the 
eclipse map, thereby yielding the most nearly axi-symmetric map that fits the data 
ES. It suppresses the azimuthal information in the default map while keeping the 
radial structure of f on scales greater than a radial blur width Ar. This seems a 
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Fig. 3 Example of the iterative fit of a data light eurve with a MEM eclipse mapping code [Bl 
Left-hand panels: data (points) and model (solid line) light curves at five different stages of the it¬ 
eration process towards convergence. Vertical dotted lines mark mid-eclipse and the ingress/egress 
of the disc centre. An horizontal dashed line depicts the uneclipsed component of the total flux. 
niter and chi2 give the cuiTent iteration number and value, respectively. Middle panels: corre¬ 
sponding eclipse maps in a logarithmic grayscale; darker regions are brighter. The notation is the 
similar to that of Fig.[^ Right-hand panels: the asymmetric component of the eclipse maps in the 
middle panels. 
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reasonable choice for accretion disc mapping because one expects the disc material 
to be roughly in Keplerian orbits, so that local departures from axi-symmetry will 
tend to be smeared away by the strong shear. This is a commonly used option and is 
also known as the default map of full azimuthal smearing. 

The full azimuthal smearing default results in rather distorted reproduction of 
asymmetric structures such as a bright spot at the disc rim. In this case, the recon¬ 
structed map exhibits a lower integrated flux in the asymmetric source region, the 
excess being redistributed as a concentric annulus about the same radial distance. By 
limiting the amount of azimuthal smearing it is possible to alleviate this effect and 
to start recovering azimuthal information in the accretion disc. Two prescriptions in 
this regard were proposed. Rutten et al. HU limited the amount of azimuthal smear¬ 
ing by averaging over a polar Gaussian weight function of constant angle Aq along 
the mapEl 


cOjk = exp 


1 ( {Rj-Rkf 

A 



( 8 ) 


while Baptista et al. mu chose to use a polar Gaussian function of constant arc 
length As through the map. 


(Ojk = exp 


1 ( {Rj-Rkf ^]k\ 


(9) 


The reconstructions in Figs.|4jb)-(d) used the default function of Eq.|^ In particular, 
the reconstructions in Figs. @1 a) and|4jb) were obtained from the same input light 
curve. This comparison illustrates that the default of limited azimuthal smearing 
provides a much better reconstruction of the central brightness distribution while 
still allowing to recover the location of a bright spot at disc rim, although with some 
azimuthal smearing. 

The concept of default map was extended with the introduction of a double de¬ 
fault function lfT9ll58lfTn . 


Dj = {DljY{D2jr ( 10 ) 

(where Dlj and D2j are separate default functions and n + ni = \), together with the 
idea of a negative default function (e.g., ni < 0). While a positive default function 
steers the MEM solution towards its intrinsic property (e.g., axi-symmetry), a nega¬ 
tive default function may be used to drive the MEM solution away from its intrinsic 
property (e.g., the presence of the undesired criss-crossed arcs) ll58l . The combina¬ 
tion of a positive axi-symmetric default function with a negative criss-crossed arcs 
default function was key to allow the recovery of accretion disc spiral structures 
with eclipse mapping ll^ fTTl . The eclipse maps shown in Figs. |2] |3 and |3 were 
obtained with this double default setup. 


^ Rj and Rf. are the distances from pixels j and k to the centre of the disc; Qji. is the azimuthal angle 
between pixels j and k\ s is the arc-length between pixels j and k. 
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Fig. 4 Effects of the default map, phase resolution (= number of data points in the light curve) and 
signal-to-noise ratio (S/N) of the input data on the eclipse mapping reconstruction. The top panel 
shows an artificial, test brightness distribution in a logarithmic grayscale. The notation is similar 
to that of Fig.[T] The middle panels show reconstructions of the test brightness distribution for 
different combination of parameters: (a) uniform default map on a light curve with N = 111 data 
points andS/Al = 100; (b) default of limited azimuthal smearing (Eq.[^,A^p = 121 and5/A^ = 100; 
(c) default of Eq.[8] Np = 121 and S/N = 20; and (d) default of Eq.[^ Np = 31 and S/N = 100. 
The lower panels show the con'esponding data (points) and model (solid line) light curves. Vertical 
dashed lines mark the ingress/egress phases of the disc centre. 


The quality of an eclipse mapping reconstruction is tied to the quality of the 
input data light curve. Specifically, the ability to recover brightness sources and 
the spatial resolution of an eclipse map depend on the phase (or, time) resolution 
and the signal-to-noise ratio {S/N) of the data light curve. Baptista & Steiner lfT4l 
provide an expression to compute the spatial resolution of the eclipse map which 
matches the phase resolution of the data light curve. Fig. |4] shows the effects of (i) 
degrading the phase resolution and (ii) lowering the S/N of the data light curve on 
the quality of the resulting eclipse map. As a reference for comparison, the eclipse 
map in Fig.|4|b) was obtained from a light curve of good phase resolution {Np=\2\ 
data points) and high signal-to-noise {S/N = 100). Fig.Uc) shows that degrading 
the S/N reduces the ability to identify fainter brightness sources. Because of the 
reduced S/N, the model light curve is not forced to follow the egress shoulder that 
signals the presence of the bright spot at disc rim. As a consequence, the eclipse 
map has only a weak asymmetry at the position of the bright spot. Fig.lUd) shows 
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that the effect of reducing the phase resolution in the light curve is to reduce the 
spatial resolution of the eclipse map, with an increase in the radial/azimuthal blur 
effects (i.e., point sources would look increasingly out of focus with decreasing 
phase resolution). The effects of low S/N seem more dramatic than those of low 
phase resolution. 


4 Error Propagation Procedures 

Because the maximum-entropy eclipse mapping is a non-linear inversion method, it 
is not possible to compute the uncertainties in the eclipse map directly from the un¬ 
certainties in the data light curve. Hence, statistical uncertainties of an eclipse map 
are usually computed with a Monte Carlo error propagation procedure ll52l . This is 
illustrated in Fig.|5] A set of N (> 20) artificial light curves is created from the orig¬ 
inal data in which the flux at each phase is varied around the true value according 
to a Gaussian distribution with standard deviation equals to the uncertainty at that 
point. The artificial light curves are fitted with the eclipse mapping code to generate 
a set of randomised eclipse maps. These are then combined to produce a “residuals” 
map by taking the pixeTto-pixel standard deviation with respect to the mean inten¬ 
sity. This yields the statistical uncertainty at each pixel. For A = 20, the uncertainty 
of the standard deviation is < 20%, sufficient to illustrate the confidence limits of 
the eclipse map, whereas N > 200 is needed in order to bring the uncertainty in 
the standard deviation down to 7% ll^ . Uncertainties obtained with this procedure 
may be used to estimate the errors in the derived radial brightness temperature and 
intensity distributions. A map of the statistical significance (or the inverse of the 
relative error) is obtained by dividing the true eclipse map by the “residuals” map 

Col EE]. 

An alternative error propagation procedure involves the use of simulations with 
the Bootstrap technique (e.g., ©)• 


5 Applications 

This section focuses on selected applications of the eclipse mapping method which 
illustrates some of the possible scientific problems that have been and may be ad¬ 
dresses with it. The reader is referred to Il34l [35] [62 [21 for more comprehensive 
reviews of the results obtained with this technique. 
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Fig. 5 Computing the statistical significance of an eclipse map from the uncertainties in the data 
light curve. Lower left panels: A set of 50 randomised light curves (dots) and comesponding eclipse 
mapping models (solid lines) superimposed in phase. The bottom panel shows these light curves 
after subtraction of the real data light curve. The scatter depicts the uncertainty in the flux at each 
phase. Top left panels: a sample of 4 of the eclipse maps from the randomised light curves, in 
a logarithmic grayscale. The notation is similar to that of Fig.|4] Top right panel: the map of the 
standard deviations with respect to the mean intensity (the “residuals” map). Lower right panel: 
the eclipse map from the real data light curve. Contour lines for S/N = 5 and 10 are overplotted; 
features in the eclipse maps are statistically significant at or above the 5 O' confidence level. 


5.1 Spectral Mapping: Spatially-resolved Disc Spectra 

The eclipse mapping method is capable of delivering spatially-resolved spectra of 
accretion discs 'when the technique is applied to time-resolved eclipse spectropho- 
tometric data, providing the best example of its ability to disentangle the light from 
different emission sources in the binary. The time-series of spectra is divided up 
into numerous 100) spectral bins and light curves are extracted for each bin. The 
light curves are then analysed to produce a series of monochromatic eclipse maps 
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Wavelength (A) 


Fig. 6 An example of spectral mapping, from HST/FOS UV-optical eclipse spectroscopy of the 
nova-like variable UXUMa. Left panel: a schematic diagram of the eclipse map, with the annular 
regions used to extract spatially resolved spectra. The disc is divided into three major azimuthal 
regions (the back side, the front side, and the quarter section containing the gas stream trajectory), 
and into a set of 6 concentric annuli with radius increasing in steps of 0.1 Ru and of width 0.05 Rn ■ 
Right panel: Spatially resolved spectra of the back region of the UX UMa accretion disc in two 
occasions (August 1994 in green; November 1994 in black). The spectra were computed for the 
set of concentric annular sections shown in the left panel. The most prominent line transitions are 
indicated by vertical dotted lines. From nni. 


covering the whole spectrum. Finally, the maps are combined to obtain the spectrum 
for any region of interest on the disc surface in. 

The spectral mapping analysis of nova-like variables ll47l |49l [TO] [T3] l29ll shows 
that their inner accretion disc is characterised by a blue continuum filled with ab¬ 
sorption bands and lines which cross over to emission with increasing disc ra¬ 
dius (Fig. |6]l. The continuum emission becomes progressively fainter and redder 
as one moves outwards, reflecting the radial temperature gradient. These high mass- 
accretion (~ 10^^ — 10^^ M© yr^^) discs seem hot and optically thick in their inner 
regions and cool and optically thin in their outer parts. The spectrum of the in¬ 
falling gas stream in UX UMa ifTOl and UU Aqr lfT3l is noticeably different from 
the disc spectrum at the same radius, suggesting the existence of gas stream “disc- 
skimming” overflow that can be seen down to ~ (0.1 —0.2) Ru¬ 
in contrast, spatially-resolved spectra of low-mass accretion (~ 10^® — 10^'^M© 
yr^*), quiescent dwarf nova show that the lines are in emission at all disc radii ll5^ . 
The observed differences between the spectra of the disc hemisphere farther from 
and closer to the LI point might be interpreted in terms of chromospheric emission 
from a disc with non-negligible opening angle (i.e., a limb-brightening effect). 

The spectrum of the uneclipsed component both for the nova-like systems and 
quiescent dwarf novae shows strong emission lines and the Balmer jump in emission 
indicating that the uneclipsed light has an important contribution from optically thin 
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gas (e.g., ini). The lines and optically thin continuum emission are most probably 
emitted in a vertically extended disc chromosphere + wind ||9l[^. This additional 
source of radiation may be responsible both for flattening the ultraviolet spectral 
slope and for Ailing in the Balmer jump of the optically thick disc spectrum, and 
might explain the historical difficulties in fitting integrated disc spectrum with stellar 
atmosphere disc models llMll^l4n[^ . 

The uneclipsed spectrum of UX UMa at long wavelengths is dominated by a 
late-type spectrum that matches the expected contribution from the mass-donor star 
ll49ll . Thus, the uneclipsed component provides an interesting way of assessing the 
spectrum of the mass-donor star in eclipsing CVs. This is a line of research yet 
to be properly explored, which could largely benefit from near- and mid-infrared 
spectroscopy of eclipses. 

Spectral mapping is also useful to isolate the spectrum of the WD, represented by 
the central pixel of the eclipse map. A stellar atmosphere model fit to the extracted 
white dwarf spectrum for the dwarf nova V2051 Oph yields the WD temperature 
and an independent estimate of the distance to the binary ll5Tl . 


5,2 Time-lapse Mapping: Dwarf Nova Outbursts 

Eclipse maps give snapshots of the accretion disc at a given time. Time-resolved 
eclipse mapping may be used to track changes in the disc structure, e.g., to assess 
variations in mass accretion rate or to follow the evolution of the surface brightness 
distribution through a dwarf nova (DN) outburst cycle. 

DN outbursts are explained in terms of either (i) the time dependent response of a 
viscous accretion disc to a burst of mass transfer from the donor star (the mass trans¬ 
fer instability model, MTIM, e.g., IflTlI ). or (ii) a limit-cycle behaviour driven by a 
thermal-viscous disc-instability (the disc-instability model, DIM, e.g., Il40ll24l '). in 
which matter progressively accumulates in a low viscosity disc (quiescence) until a 
critical surface density is reached at a given radius, causing a heating wave to switch 
the disc to a high viscosity regime (outburst) that allows the gas to diffuse rapidly 
inwards and onto the white dwarf. There is a set of distinct predictions from these 
models that might be critically tested with time-lapse eclipse mapping 0. 

An example of time-lapse eclipse mapping is given in Fig.|7] Eclipse maps cov¬ 
ering the full outburst cycle of the dwarf nova EX Dra Q show the formation of a 
one-armed spiral structure at the early stages of the outburst and reveal how the disc 
expands during the rise until it fills most of the primary Roche lobe at maximum 
light. During the decline stage, the disc becomes progressively fainter until only a 
small bright region around the WD is left at minimum light. The evolution of the 
radial temperature distribution shows the presence of an outward-moving heating 
wave during rise and of an inward-moving cooling wave in the decline. 

The temperatures of the EX Dra outbursting disc are above those expected for 
the hot, ionised outbursting disc gas in the DIM framework ll^l40ll . and the cooling 
wave decelerates as it travels inwards Q. Since these results are consistent with both 
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binary phase 


log( radius [R/RLl] ) 


Fig. 7 Time-lapse eclipse mapping of the dwarf nova EX Dra along its outburst cycle. Left-hand 
panels: Data (dots) and model (solid line) light curves in (a) quiescent, (b) early rise, (c) late rise, 
(d) outburst maximum, (e) early decline, and (f) late decline stages. Vertical dotted lines mark 
ingress/egress phases of disc centre. Middle panels: eclipse maps in a false colour blackbody loga¬ 
rithmic grayscale. The notation is similar to that of Fig.|3 The numbers in parenthesis indicate the 
time (in days) elapsed since outburst onset. Right-hand panels: azimuthal-averaged radial bright¬ 
ness temperature distributions for the eclipse maps in the middle panels. Dashed lines show the l-a 
limit on the average temperature for a given radius. A dotted vertical line depicts the radial position 
of the bright spot in quiescence; vertical ticks mark the position of the outer edge of the disc (red) 
and the radial position at which the disc temperature falls below 11000 K (blue). Steady-state disc 
models for mass accretion rates of log M= —8.0 and —9.0 Mq yr“* are plotted as dotted lines for 
comparison. Numbers in parenthesis list the integrated disc luminosity. Erom (7|. 
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DIM and MTIM, they offer no power to discriminate between these models. How¬ 
ever, the same eclipse mapping study reveals that the radial temperature distribution 
in quiescence follows the T law of opaque steady-state discs (Fig.|3a)), 

and that this quiescent state is reached only 2 d after the end of the outburst. Both 
results indicate that the viscosity of its quiescent disc is as large as in outburst - in 
contradiction with DIM and in agreement with predictions of the MTIM. In addi¬ 
tion, the early rise map shows evidence of enhanced gas stream emission, indicative 
of an enhanced mass transfer rate at this early outburst stage (Fig.|2lb)). The inte¬ 
grated disc luminosity at early rise (1.4±0.31,0) is comparable to that in quiescence 
(1.2 ± 0.1 Lq), and is not enough to support the idea that the enhanced mass transfer 
could be triggered by an increased irradiation of the mass donor by the accretion 
disc. This led to the conclusion that the observed enhanced mass transfer at early 
rise is not a consequence of the ongoing outburst, but its cause - suggesting that the 
outbursts of EX Dra are powered by bursts of mass transfer 0. 

Time-lapse eclipse mapping studies were also performed for the dwarf novae 
ZCha ia, OYCar lE]El, IP Peg ([T9], V2051 Oph US, and HTCas ll^l26l. 


5.3 Flickering Mapping: Revealing the Disc Viscosity 


Flickering is the intrinsic brightness fluctuation seen in light curves of T Tau stars 
1311, mass-exchanging binaries ||T||22] and active galactic nuclei lIZTl . Optical stud¬ 
ies suggest there are two different sources of flickering in CVs: (i) the stream-disc 
impact region (possibly because of unsteady mass inflow or post-shock turbulence 
ll^ l55l l and/or (ii) turbulent inner disc regions (possibly because of magneto¬ 
hydrodynamic (MHD) turbulence, unsteady WD accretion or events of magnetic re¬ 
connection at the disc atmosphere ll28ll2n[J7l ). With its ability to spatially-resolve 
and to disentangle different sources, flickering mapping has been a useful tool to 
confirm and extend this scenario. As an added bonus, if the disc-related flickering 
is caused by MHD turbulence, one may infer the radial run of the disc viscosity 
parameter ass El from the relative flickering amplitude, (Gd/D) ll28ll . 


a, s{R) ^0.23 


R 

(yoiR) 

50H 

0.05D{R) 


( 11 ) 


where H is the disc scale height. 

From a large, uniform ensemble of light curves of a given CV it is possible to 
separate the steady-light component (the average flux in a given phase bin), low- and 
high-frequency flickering amplitudes (the scatter with respect to the average flux, 
computed with the ensemble lITSIl and single ll2Tl methods) as a function of binary 
phase, to derive corresponding maps of surface brightness distributions from their 
eclipse shapes and, thereafter, to compute the radial run of the relative amplitude of 
the disc flickering component liSil^. 
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Flickering mapping of the DN V2051 Oph reveals that the low-frequency flick¬ 
ering arises mainly in an overflowing gas stream and is associated with the mass 
transfer process. The high-frequency flickering originates in the accretion disc and 
has a relative amplitude of a few percent, independent of disc radius and brightness 
level, leading to large values ttss — (0.1 — 0.2) at all disc radii ||5l . 

Figure[8]shows the results of the eclipse mapping analysis of an extensive data set 
of optical light curves of the dwarf nova HT Cas 0. These observations frame a 2 d 
transition from a low state (largely reduced mass transfer rate) back to quiescence, 
allowing the application of both time-lapse and flickering mapping techniques to 
derive independent estimates of ttss- In the low state, the gas stream hits the disc at 
the circularisation radius Rare ISQIj the accretion disc has its smallest possible 
size. The disc fast viscous response to the onset of mass transfer, increasing its 
brightness and expanding its outer radius at a speed v ~ +0Akms^^, implies ass — 
0.3 — 0.5. The newly added disc gas reaches the WD at disc centre soon after mass 
transfer recovery (~ 1 d), also implying a large disc viscosity parameter, a^s — 0.5. 
Flickering mapping reveal a minor, low-frequency BS-stream flickering component 
in the outer disc, plus a main disc flickering component the amplitude of which rises 
sharply towards disc centre (Fig.|8]l, leading to a radial dependency ass{R) R^^ 
with ttss >0.1 for R < Rare - in agreement with the time-lapse results. 

A similar analysis was performed for the DN V4140 Sgr a. Eclipse mapping in 
quiescence indicate that the steady-light is dominated by emission from an extended 
disc with negligible contribution from the WD, suggesting that efficient accretion 
through a high-viscosity disc is taking place. Flickering maps show an asymmet¬ 
ric source at disc rim (BS-stream flickering) and an extended central source (disc 
flickering) several times larger in radius than the WD. Unless the thin disc approxi¬ 
mation breaks down, the relative amplitude of the disc flickering leads to large ttss’s 
in the inner disc regions (~ 0.2 — 1.0), which decrease with increasing radius. 

In contrast, in the nova-like UU Aqr optical flickering arises mainly in tidally- 
induced spiral shocks in its outer disc 0. Assuming that the turbulent disc model 
applies, its disc viscosity parameter increases outwards and reaches a^j ~ 0.5 at the 
position of the shocks, suggesting that they might be an effective source of angular 
momentum removal of disc gas. 


5.4 3D Eclipse Mapping: Disc Opening angle & Superhumps 

Standard eclipse mapping assumes a simple flat, geometrically thin disc model. Real 
discs may however violate this assumption in the limit of high M. Disc half-opening 
angles of )3 > 4" are predicted forM> 5 x lO^^MQjr^^ ll43l . At large inclinations 
(i > 80'’) this may lead to artificial front-back asymmetries in the eclipse map be¬ 
cause of the different effective areas of surface elements in the near and far sides 
of a flared disc as seen by an observer on Earth. Eurthermore, the assumption that 
the emitted radiation is independent of orbital phase implies that any out-of-eclipse 
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Fig. 8 Left-hand panels: data (green dots) and eclipse mapping model (black lines) disc light 
curves, and the extracted light curve of the WD (red lines). Top and middle panels show, respec¬ 
tively, the low state and quiescence steady-light component, while the lower panel shows the high- 
frequency (/ > 3.3mHz) flickering curve. Middle panels: the corresponding eclipse maps in a 
logarithmic grayscale. The notations is similar to that of Fig.[T] Dashed circles mark the outer disc 
radius in each case. Right-hand panels: azimuthal-averaged radial intensity (top, middle) and flick¬ 
ering amplitude in percent (bottom) distributions for the eclipse maps in the middle panels. Dashed 
lines show the 1-CT limit on the average value for a given radius. 


brightness change (e.g., orbital modulation due to BS anisotropic emission) has to 
be removed before the light curve can be analysed!!. 

A step to overcome these limitations is to go beyond the standard assumptions 
allowing the eclipse mapping surface to become three-dimensional. This leads to 
3D eclipse mapping m- For example, with the inclusion of a disc rim in the 
eclipse mapping method, the out-of-eclipse modulation can be modelled as the fore¬ 
shortening of an azimuthal-dependent brightness distribution in the disc rim im. 
This procedure allows one to recover the azimuthal (phase) dependency of the BS 
emission. Moreover, motivated by the front-back asymmetry that appeared in the 
flat-disc map and by the difficulties in removing the asymmetry with the assumption 
of an uneclipsed component, na introduced a flared disc in their eclipse mapping 


^ Thi.s is usually done by fitting a spline function to the phases outside eclipse, dividing the light 
curve by the fitted spline, and scaling the result to the spline function value at phase zero (e.g., |^). 
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of ultraviolet light curves of Z Cha at outburst. They found that the asymmetry van¬ 
ishes and the disc is mostly axi-symmetric for a disc half-opening angle of j3 = 6". 

In 3D eclipse mapping, the mapping surface usually consists of a grid of A^l 
pixels on a conical surface centred at the WD position and inclined at a half-opening 
angle j3 with respect to the orbital plane, plus a circular rim of N2 pixels orthogonal 
to the orbital plane at a distance Rd (< Rli) from disc centre. An entropy landscape 
technique ll4^ may be used to cope with the extra degree of freedom that comes 
along with the additional geometry parameters p and Rd- Simulations show that 
if an eclipse mapping reconstruction is performed with the wrong choice of j3 and 
Rd, the code develops artefacts in the brightness distribution in order to compensate 
for the incorrect parameters. Eclipse maps with these spurious, additional structures 
have lower entropy than the map with the correct choice of j3 and Rd- Therefore, the 
best-ht geometry can be found by searching the space of parameters for the pair of 
{P,Rd) values of highest entropy (Fig.|9] left-hand panels). 

An example of 3D eclipse mapping is shown in Fig.|9] where the technique was 
applied to study the evolution of the accretion disc surface brightness of the DN 
V4140Sgr in a superoutburst 0. The entropy landscape analysis indicates that the 
accretion disc is geometrically thin both in outburst (j3 = 1.0“ ± 0.5°) and in quies¬ 
cence (P ~ 0.5° ± 0.5"); it hlls the primary Roche lobe in outburst and progressively 
shrinks to about half this size in quiescence. They also hnd that the disc is elliptical 
in outburst and decline, with an eccentricity e = 0.13. At both outburst stages, the 
disc orientation is such that superhump maximum occurs when the mass-donor star 
is aligned with the bulge of the elliptical disc. This lends observational support for 
the tidal resonance instability model of superhumps ll^l^l42l . 

Additional observational support for the presence of elliptical precessing discs 
in CVs come from the eclipse mapping analyses of light curves of the DN ZCha in 
superoutburst ll44l and of the permanent superhumper V348 Pup ll46l . 


6 Summary 

Eclipse mapping is a unique, powerful technique to investigate: 

• Accretion disc spectra (allowing one to separate the disc atmosphere emission at 
different distances from disc centre, and to isolate the spectrum of the WD, BS, 
a possibly outflowing disc wind, and even the faint, red mass-donor star); 

• Accretion disc structures (such as gas stream overflow, tidally-induced spiral 
shocks, elliptical precessing discs, and magnetic accretion curtains); 

• Time evolution of accretion discs (tracing the mass and angular momentum re¬ 
distribution throughout dwarf nova outbursts, or to follow brightness changes 
caused by lighthouse effects in discs with fast spinning, magnetic WDs); 

• Accretion disc viscosity (either via flickering mapping or by measuring the ve¬ 
locity of transition waves during dwarf nova outbursts with time-lapse mapping). 
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Fig. 9 3D eclipse mapping of V4140Sgr. Left panels: entropy landscape for the outburst (top), 
decline (middle) and quiescence (bottom) light curves in a blackbody false-colour scale. Re¬ 
gions in yellow have higher entropy. A red cross marks the combination of disc half-opening 
angle j8 and rim radius Rj of highest entropy. Middle panels: data (dots) and model (solid line) 
light curves for the pair of (/3,Rrf) values of highest entropy in each case. Vertical dotted lines 
mark the ingress/egress phases of the WD, while vertical tick marks depict the phases of eclipse 
ingress/egress. Right panels: corresponding eclipse maps in a logarithmic grayscale. The notation 
is similar to that of Fig.[2] Dashed circles show the best-fit disc rim radius. From a 
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